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I. Introduction 



In the early 1970's Bekenstein argued that black holes have entropy [1] 
and Hawking showed that black holes have temperature [2] . A few years later 
these conclusions were apparently confirmed by Gibbons and Hawking [3] 
who used path integral methods to evaluate the black hole partition function. 
However, the original analysis of Gibbons and Hawking is incorrect. The 
difficulty arises because, in the language of ordinary thermodynamics, black 
holes are unstable and have a negative heat capacity. On the other hand, 
any system that can be described by a partiton function is necessarily stable 
and has a positive heat capacity. 

The path integral approach of Gibbons and Hawking is also somewhat 
mysterious regarding the origin of black hole entropy. A number of authors 
have speculated that the large entropy of a black hole should be associated 
with a large number of internal states, hidden by the horizon, that are con- 
sistent with the few external parameters that characterize the black hole. 
Others have argued that black hole entropy can be associated with a large 
number of possible initial states that can collapse to form a given black hole. 
The path integral analysis does not support either of these views in an obvi- 
ous way. The key element in the path integral calculation is the action of a 
static Euclidean black hole geometry. A static Euclidean black hole has no 
horizon, no interior region, and no information about the collapse process. 
The path integral calculation of the partition function seems to hint at a 
topological explanation for black hole entropy. 

The first objective of this article is to show that the black hole partition 
function can be placed on a firm logical foundation by enclosing the black hole 
in a spatially finite "box" or boundary. The presence of the box has the effect 
of stabilizing the black hole and yields a system with a positive heat capacity. 
Without a finite box, or some other mechanism for stabilizing the black hole 
(such as a negative cosmological constant [4]), the partition function does 
not exist. It is, nevertheless, common practice to fix boundary conditions at 
infinity and treat the zero-loop approximation to the path integral as if it 
were a real partition function. Such an identification is logically unfounded 
and, at any rate, cannot be extended beyond the zero-loop approximation. 

The second objective of this article is to explore the origin of black hole 
entropy. This is accomplished through the construction of a path integral 
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expression for the density matrix for the gravitational field, and through an 
analysis of the connection between the density matrix and the black hole 
density of states. (The density matrix for a black hole has been studied 
previously in Ref. [5].) Our results suggest that black hole entropy can be 
associated with an absence of certain "inner boundary information" for the 
system. 

We begin in Sec. 2 with a review of the motivations for enclosing the 
black hole in a box. In Sec. 3 we analyze the action for the gravitational 
field in the presence of a spatially finite boundary. Section 4 contains a 
discussion of the relationship between Lorentzian and Euclidean notation for 
the gravitational action. Formal path integral expressions for the density 
matrix, density of states, and partition function are constructed in Sec. 5. In 
Sec. 6 we use these results to calculate the density of states for a black hole 
and we discuss the origin of black hole entropy. 

II. Finite Boundaries 

Hawking's analysis [2] shows that the temperature of a black hole, as 
measured at spatial infinity, equals the surface gravity divided by 2ir. For a 
Schwarzschild black hole of mass M, it follows that the inverse temperature 
[3 at infinity is 8nM. On the other hand, the standard thermodynamical 
definition of inverse temperature [6] is f3 — dS(E)/dE, where S(E) is the 
entropy function and E is the thermodynamical internal energy. If the mass 
at infinity M and the internal energy E are identified, then the relationship 
dS(E)/dE = 8nM can be integrated to yield S(E) = AttE 2 (plus an additive 
constant). This result is in complete agreement with the prediction made by 
Bekenstein [1] that a black hole has entropy proportional to the area of its 
event horizon. 

The black hole entropy S(E) = AttE 2 is a convex function of E. This 
is characteristic of an unstable thermodynamical system [6]. In this case, 
the instability arises because energy and temperature are inversely related 
for black holes. Thus, if fluctuations cause a black hole to absorb an extra 
amount of thermal radiation from its environment, its mass will increase and 
its temperature will decrease. The tendency then is for the cooler black hole 
to absorb even more radiation from its hotter environment, causing the black 
hole to grow without bound. 
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These results can be reformulated within the context of statistical me- 
chanics. First consider the canonical partition function Z(f3) for an arbitrary 
system. In general Z((3) is a sum over quantum states weighted by the Boltz- 
mann factor e~ l3E . If v(E) is the density of quantum states with energy E, 
then 

Z((3) = j dEu(E)e^ E . (2.1) 

The partition function can also be expressed as 

Z((3) = J dEe- I{E) , (2.2) 

where the "action" is defined by 1(E) = f3E — S(E) and the entropy function 
S(E) is the logarithm of the density of states: S(E) = \nu(E). (Dimension- 
ful constant factors can be included in the logarithm as necessary.) The 
integral over E can be evaluated in a steepest descents approximation by 
expanding the action 1(E) to quadratic order around the stationary points 
E*(j3), which satisfy 
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The Gaussian integral associated with a stationary point E* will converge if 
the second derivative of the action at E* is positive: 
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This condition shows that the entropy S(E) should be a concave function at 
the extremum E* in order for the Gaussian integral to converge. 

A further significance of the condition (2.4) can be seen as follows. In the 
steepest descents approximation, the expectation value of energy is (E) = 
-dlnZ/dp « E* and the heat capacity is C = d(E)/dj3~ 1 « dE*/dp~\ By 
differentiating Eq. (2.3) with respect to (3, we find 
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Thus, we see that in the steepest descents approximation the convergence of 
the integral for the canonical partition function is equivalent to the thermo- 
dynamical stability of the system (the concavity of the entropy), which in 
turn is equivalent to the positivity of the heat capacity. 

For the black hole in particular, the entropy S(E) = AttE 2 is not a 
concave function of the internal energy and the integral for the partition 
function diverges. This can be seen in the path integral formalism as well. 
In that case, the partition function is expressed as a functional integral over 
Euclidean metrics with periodic time, where the time period equals the in- 
verse temperature at infinity [3]. The action is extremized by a Euclidean 
black hole metric, but the integration about that stationary point diverges. 
Formally, such a divergent Gaussian integral yields an imaginary result. The 
Euclidean black hole extremum is properly interpreted as an instanton that 
dominates the semiclassical evaluation of the rate of black hole nucleation 
from flat space at finite temperature [7]. 

The preceeding discussion shows that the canonical partition function 
Z{($) characterizes the thermal properties of thermodynamically stable sys- 
tems. For unstable systems Z(/3) can give information concerning the rate 
of decay from a quasi-stable configuration (such as "hot flat space" in the 
black hole example), but it cannot be used to define thermodynamical prop- 
erties such as expectation values, fluctuations, response functions, etc. Thus, 
before the partition function can be used as a probe of black hole thermo- 
dynamics, it is first necessary to stabilize the black hole. It was recognized 
in Ref. [8] that a black hole is rendered thermodynamically stable by enclos- 
ing it in a spatially finite "box" or boundary whose walls are maintained at 
a finite temperature. In this case the energy and the temperature at the 
boundary are not inversely related because of the blueshift effect of tempera- 
ture in a stationary gravitational field [9] . If the black hole absorbs an extra 
amount of thermal radiation from its environment, its energy will increase 
thereby increasing the gravitational blueshift. Although the temperature of 
the black hole as measured at infinity is decreased by such a fluctuation, the 
temperature as measured at the boundary can increase due to the enhanced 
blueshift effect. The hotter black hole will then give up its excess energy and 
return to its stable equilibrium configuration. 

The stabilizing effect of a finite box can be confirmed by the following 
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simple analysis. Consider a Schwarschild black hole of mass M surrounded 
by a spherical boundary of radius R. The inverse temperature at infin- 
ity is 8nM, so the inverse temperature at the boundary is blueshifted to 
(3 = 8nM^Jl — 2M/R [9]. On the other hand, the inverse temperature is 
defined by (3 = dS(E)/dE, where again S(E) is the entropy as a function 
of internal energy E. Now, the entropy of the black hole, at least in the 
zero-loop ("classical") approximation, depends only on the black hole size 
and is unaffected by the presence, absence, or proximity [8] of a finite box. 
Thus, we have S(E) = AttM 2 as before. By equating the two expressions for 
inverse temperature we find 

SnMJ^m^^fl. (2.7) 

oE 

In this case the energy E and the mass M as measured at infinity do not 
coincide. Equation (2.7) can be integrated to yield [8,10] 

E = R- Ry/l - 2M/R , (2.8) 

where, for convenience, the integration constant has been chosen so that 
E — > M in the limit R — > oo with M fixed. The significance of this expression 
can be seen by expanding E in powers of GM/R (where Newton's constant 

G is set to unity), with the result E = M + M 2 /(2R) H . This shows 

that the internal energy inside the box equals the energy at infinity M minus 
the binding energy —M 2 /(2R) of a shell of mass M and radius R. The 
binding energy —M 2 /(2R) is the energy associated with the gravitational 
field outside the box [10]. Also observe that the internal energy takes values 
in the range < E < R. 

By solving Eq. (2.8) for M as a function of E, we obtain the entropy 
function 

S(E)=AnE 2 (l-E/(2R)) 2 . (2.9) 

First, note that the derivative dS / dE is a concave function of E that vanishes 
at the extreme values E = and E = R. It follows that dS/dE has a 
maximum (3 cr . For (3 > f3 cr the equation (3 = dS/dE has no solutions for 
E. On the other hand, for f3 < (3 cr , there are two solutions E\ and E 2 . At 
the larger of these two solutions, say, E 2 , the second derivative d 2 S/dE 2 is 
negative and the stability criterion (2.4) is satisfied. At the smaller of these 
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two solutions, Ei, the second derivative d 2 S/dE 2 is positive and the stability 
criterion (2.4) is violated. These considerations indicate that for a small box 
at low temperature {(5 > fl cr ), the equilibrium configuration consists of flat 
space. For a large box at high temperature {(5 < f3 cr ), the stable equilibrium 
configuration consists of a large black hole with energy E 2 . The unstable 
black hole with energy E\ is an instanton that governs the nucleation of 
black holes from flat space. In the limit R — > oo, the stable black hole 
configuration is lost and only the instanton solution survives [8,4]. 

In the following sections, we will develop the formal functional integral 
expressions for the density of states, the density matrix, and the canonical 
partition function for the gravitational field with a spatially finite bound- 
ary B. For the partition function in particular, the inverse temperature will 
be fixed as a boundary condition on B. The use of a finite boundary has 
the following very important consequence: Because gravitational fields cause 
temperature to redshift and blueshift, one must allow for the temperature 
to be fixed to different values at different points on B. In other words, 
gravitating systems in thermal equilibrium are not characterized by a single 
temperature but instead by a temperature field on the boundary of the sys- 
tem [11,12]. Correspondingly, the partition function is actually a. functional 
Z[f3] of the inverse temperature field on B. For some problems, such as the 
Schwarzschild black hole in a spherical box discussed above, it is possible 
to choose the temperature to be a constant on B. In those cases B coin- 
cides with an isothermal surface for the system. However, experience with 
the Kerr black hole shows that this must be viewed as a particular choice 
of boundary conditions, not the most general choice. What happens in the 
Kerr case [11] is that the angular velocity of the black hole with respect to 
observers who are at rest in the stationary time slices enters as a "chemical" 
potential conjugate to angular momentum. It turns out that the constant 
temperature surfaces and the constant angular velocity surfaces do not co- 
incide. Therefore it is necessary to allow for some thermodynamical data, 
either the temperature or the chemical potential or both, to vary across the 
boundary. This conclusion might seem disturbing at first, since traditionally 
one of the purposes of thermodynamics has been to provide a characteriza- 
tion of systems in terms of only a few parameters. Note, however, that the 
thermodynamic parameters such as j3 are only given on a two-surface B. 
No assertion is made about the values of (3 inside B. Thus the generalized 



6 



thermodynamical description required is still far simpler than that of the 
full configuration of the system. The thermodynamical formalism that re- 
sults from a generalization to non-constant thermodynamical data also has a 
number of compelling features. In particular, the thermodynamical data are 
brought into direct correspondence with the canonical boundary data, and in 
the process an intimate connection between thermodynamics and dynamics 
is revealed [11-13]. 

III. The Action 

Assume that the spacetime manifold M. is topologically the product of 
a spacelike hypersurface and a real line interval, Ex J. The boundary of £ 
is denoted <9£ = B. The spacetime metric is with associated curvature 
tensor 1Z^ vap and derivative operator V M . The boundary of M, dM, con- 
sists of initial and final spacelike hypersurfaces t' and t", respectively, and a 
timelike hypersurface 3 B = B x I joining these. The induced metric on the 
spacelike hypersurfaces t' and t" is denoted by hij, and the induced metric 
on 3 B is denoted by 7^. 2 

Consider the gravitational action 
S 1 = — [ d A x^g~(1l - 2A) + - f d 3 xVhK - - [ d 3 x^Q . (3.1) 

2k Jm K Jt> K J' A B 

Here, k is 8tt times Newton's constant and A is the cosmo logical constant. 
For simplicity, we have omitted matter contributions to the action. The 
symbol f t , d 3 x denotes an integral over the boundary element t" minus an 
integral over the boundary element t'. The function K is the trace of the 
extrinsic curvature Kij for the boundary elements i! and t", defined with 
respect to the future pointing unit normal u^. Likewise, G is the trace of the 
extrinsic curvature G^- of the boundary element 3 B, defined with respect to 
the outward pointing unit normal n^. 

Under variations of the metric the action (3.1) varies according to 

SS 1 = (terms that vanish when the equations of motion hold) 

2 We use latin letters i, j, k, ... as indices both for tensors on 3 B and for tensors on a 
generic hypersurface S. The two uses of such indices can be distinguished by the context 
in which they occur. 
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+ 



ft" r | fB" 

J d 3 xP l] 5h ij + J d 3 xTr t3 6jij- -J d 2 x^5a . (3.2) 



The coefficient of Shij in the boundary terms at t' and t" is the gravitational 
momentum 

P« = —Vh(Kh ij - K ij ) . (3.3) 
2/t 

Likewise, the coefficient of 5^ in the boundary term at 3 B is 

ir ij = -^-V^f(Qf j - Q ij ) ■ (3.4) 

Equation (3.2) also includes integrals over the "corners" B" = t" fl 3 B and 
B' — t' n 3 I? whose integrands are proportional to the variation of the "angle" 
a = sinh _1 (-u • n) between the unit normals of the hypersurfaces t" and t' 
and the unit normal of 3 B [14,15]. The determinant of the two-metric on 
B' or B" is denoted by a. 

The action S* 1 yields the classical equations of motion when the induced 
metric on 3 B, and t" and the angle a at B' and B" are held fixed in the 
variational principle. In general, the functional S = S 1 — S°, where S° is 
a functional of the fixed boundary data, also yields the classical equations 
of motion. For simplicity we define S° to be a functional of 7^ only. The 
variation 5S then differs from 5S 1 of Eq. (3.2) only in that 7r lJ is replaced by 

Now foliate the boundary element 3 B into two-dimensional surfaces B 
with induced two-metrics a a b- The three-metric 7^ can be written according 
to the familiar Arnowitt-Deser-Misner decomposition as 

7y dx i dx j = -N 2 dt 2 + a ab (dx a + V a dt) (dx b + V b dt) , (3.5) 

where N is the lapse function and V a is the shift vector. The corresponding 
variation of 7^ is [10] 

5 ll3 = (-2u iUj /N)SN + (-2a a{lUj) /N)5V a + {ofc^Sa* , (3.6) 

where Ui is the unit normal of the slices B and o % a = 5 l a projects covariant 
tensors from 3 B to the slices B. With this result, the contribution to the 
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variation of S from the boundary element 3 B becomes 



SS\ 3B = I d 3 x{^ - (6S°/6>y ij ))6>y ij 

= J d 3 x^(-e5N + 3a 5V a + (N/2)s ab 5a ab ) , (3.7) 

where the coefficients of the varied fields are defined by 

2 1 SS° 

2 1 5S° 

2 2^9° 

SClb = TT-r^^i- AT I — JT ■ ( 3 - 10 ) 

The leading terms in Eqs. (3.8)-(3.10) can be rewritten in terms of an extrin- 
sic curvature k ab that is defined by parallel transporting the unit normal 
of 3 B across a two-dimensional slice B. Thus, k ab is the extrinsic curvature of 
B considered as the boundary B = <9S of a spacelike hypersurface £ whose 
unit normal w M is chosen orthogonal to n^. (The case in which u^n^ ^ 
has been considered in Refs. [16,17], and will be treated in Ref. [15].) Also 
let P lJ denote the gravitational momentum for the hypersurfaces E that are 
"orthogonal" to 3 B, and let = u u V u u^ denote the acceleration of the unit 
normal for this family of hypersurfaces. The resulting expressions are [10] 

e = —k — e , (3-11) 

K 

Ji = -^=<r ij P jk n k -(j )i, (3-12) 
s ab = -(k ab + (nX - k)a ab ) - (s ) ab . (3.13) 

K 

In these equations, we have expressed the terms proportional to the functional 
derivatives of S° as e , (j )i, and (s {) ) ab . We will assume that S* is a linear 
functional of the lapse N and shift V a , so that e and (j„); are functionals of 
the two-metric a ab only [10]. Also note that the indices in Eq. (3.12) refer 
to the hypersurface S. Thus, ji = j a cr? where a l a = 5 l a projects tensors from 
£ to B, and = crfa a j. 
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From its definition through Eq. (3.7), -Joe equals minus the time rate 
of change of the action, where changes in time are controlled by the lapse 
function N on 3 B. Thus, e is identified as an energy surface density for 
the system [10]. Likewise, we identify j$ as the momentum surface density 
and s ab as the spatial stress [10]. The fact that the relevant stress-energy- 
momentum quantities live on the boundary of space is a consequence of 
the fact that the gravitational field on the boundary B "knows" about the 
entire matter-energy content of space interior to B. Also note that the total 
energy E can be defined by integrating \fae over the two-surface B. The 
result is precisely the energy discussed in Sec. 2. In particular, if B is a 
spherical boundary in a static time slice of a Schwarzschild black hole, the 
total energy E is given by Eq. (2.8). (This assumes S* is chosen such that 

E = -1/(47^).) 



IV. Lorentzian and Euclidean Notation 

The action S discussed in the previous section can be viewed as a func- 
tional of complex metrics. In particular S is the action for both Lorentzian 
and Euclidean metrics. We adopt the point of view that, strictly speaking, 
there is no distinction between the "Lorentzian action" and the "Euclidean 
action" , or between the "Lorentzian equations of motion" and the "Euclidean 
equations of motion" . Of course, a particular solution of the classical equa- 
tions of motion can be Lorentzian or Euclidean. But for the action functional 
itself, the only distinction between Lorentzian and Euclidean is simply one 
of notation. In Sec. 3 we have used what might be called Lorentzian nota- 
tion: the action S is defined with the convention that exp(iS) is the phase 
in the path integral; the volume elements for Ai and 3 B are written as a/— g 
and respectively; the lapse function associated with the foliation of M. 

into hypersurfaces E is defined by N = ^ — l/g tt . (The lapse function that 
appears in Eq. (3.5) is the restriction of this spacetime lapse to the bound- 
ary element 3 B. It is defined by N = ^J—l/j tt .) Therefore S, y/^^g, \J— 7, 
and N are real for Lorentzian metrics and imaginary for Euclidean metrics. 
We can re-express the action in Euclidean notation by making the following 
changes. Define a new action functional by I[g] = —iS[g] so that the phase 
in the path integral is given by exp(— I). Also rewrite the volume elements 
for M and 3 B as yjg = i^J^g and y/^y = iyj—'y, respectively, and define a 
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new lapse function by N = J l/g u = iy—l/g tt = iN.f\ 

It is also convenient to redefine the timelike unit normal of the slices E. 
In Lorentzian notation, the unit normal is defined by = —NSj^ and satisfies 
u • u — — 1. A new unit normal is defined by u M = N5* = iNSt = —iu^, and 
satisfies u-u — +1. In some contexts it is also useful to define a new extrinsic 
curvature K^ u in terms of the normal u^. K^ u is related to the old extrinsic 
curvature by = -(6° - u a u^)V a u v = i(6° + ifu^V ' c u v = -iK^ v . 
In turn, K^ v can be used to define a new gravitational momentum P*- 7 that 
is related to the momentum of Eq. (3.3) by P % i = —iP l K We will, however, 
continue to use the old notations Kij and P % K 

In addition to the notational changes described above, we will also define 
a new shift vector by Vi = igu = %V%- This notation is a departure from the 
standard Euclidean notation in the sense that Vi is imaginary for Euclidean 
metrics. One of the motivations for this change is the following. Apart from 
surface terms, the gravitational Hamiltonian is a linear combination of con- 
straints built from the gravitational canonical data with the lapse function 
and shift vector as coefficients. In conjunction with the new notation N, 
Vi for the lapse and shift, we choose to continue to denote the gravitational 
canonical data by hij, P % '\ as mentioned above. Then the constraints are 
unaffected by the change in notation, and the Hamiltonian can be written 
as H[N,V] = —iH[N,V]. The overall factor of —i that appears in this 
relationship is precisely what is required for the connection between the evo- 
lution operator [e~ %m in particle mechanics) and the density operator (e~ H/S 
in ordinary statistical mechanics). When the gravitational field is coupled 
to other gauge fields, such as Yang-Mills or electrodynamics, it is natural 
to redefine the Lagrange multipliers associated with the gauge constraints as 
well [11]. 

With our new notation, Eq. (3.1) becomes 
I l = -—( d 4 aVHft - 2A) - - /* d 3 xVhK+- [ a?xJ^Q (4.1) 

2K Jm V ' K JH K i 3 B V 

3 A bit of care is required in defining the square roots. For example, the appropriate 
definition of \/—g is obtained by taking the branch cut in the upper half complex plane, 
say, along the positive imaginary axis. Then the imaginary part of \f— g is negative. 
Correspondingly, the appropriate definition of ^fg is obtained by taking the branch cut 
along the negative imaginary axis. Then the imaginary part of ^fg is positive. 
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and Eq. (3.2) becomes 

5I 1 = (terms that vanish when the equations of motion hold) 

rt" r 1 pB" 

— % J d 3 xP l] 5h ij -i J 3B d3x7ll36 ^ + -j B/ d 2 x^5a . (4.2) 

Here, we have defined a = cos -1 (iZ-n) so that 5a = iSa. Thus, a is the angle 
between the unit normals u and n of the boundary elements t" (or t') and 
3 B. The full action / = —iS differs from I 1 by a term 1° = — iS° that is a 
functional of the metric 7jj on 3 B. The contribution to 5/ from the boundary 
element 3 B is obtained from Eq. (3.7). Putting this together with Eq. (4.2) 
yields 

51 = (terms that vanish when the equations of motion hold) 

rt" 1 rB" 

— i I d 3 x P %3 6hij + - d 2 x^/a5c 

Jt' K JB> 

+ J 3B d 3 xV^ (eSN - j a 5V a - (N/2)s ab 5a ab ) . (4.3) 

The stress-energy-momentum quantities e, ji, and s ab from Eqs. (3.11)- 
(3.13) are given by 

< = lk -^m- <*■*) 

K -JdSN 

]i = ~7f" P nt + 7Ssv^ (4 ' 5) 

In Eq. (4.6) we have written the n M component of the acceleration of the unit 
normal as n^a^ = (n l diN)/N = (n l diN)/N. Also note that the action / 
in canonical form is [10,12] 

/ = j^d A x (-iP+'hij + NH + \T7ii) + d 3 xy/a (Re - V a j a ) . (4.7) 

Here, 7i and Tii are the Hamiltonian and momentum constraints. 

Since 1° (like S°) is functionally linear in the lapse N and shift V a , then e 
and ji are functions of the canonical variables only. Moreover, s ab depends on 



12 



the lapse and shift only through ratios such as (diN)/N and V a /N. These 
observations are important for the following reason. Consider the metric 
g L for a stationary Lorentzian geometry, written in stationary coordinates. 
Given such a metric, we can compute the stress-energy-momentum quanti- 
ties e, ji, and s ab for the surface 3 B that is defined by taking a two-surface 
B within one of the stationary time slices and extending this surface along 
the stationary time direction. Now consider the generally complex metric 
g c obtained by substituting t — > —it in the stationary Lorentzian metric 
g L . Under this transformation the lapse function N and shift vector V 1 are 
changed from real to imaginary, and the canonical variables hij, P 1 ^ are in- 
variant. That is, the substitution t — > —it generates a complex metric g c 
with N, V" 1 , hij, and P* J real. (The metric will be Euclidean if V % = 0.) 
If we now compute e, ji, and s ab for the surface 3 B for the complex metric, 
the results will be the same as the results for the original Lorentzian metric. 
This observation is the key to understanding the relationship between the 
results of a functional integral over complex (or Euclidean) metrics and the 
real Lorentzian spacetime that it describes [12]. For example, suppose we 
compute the partition function as a functional integral and find that it is 
extremized by a certain stationary complex (or Euclidean) black hole metric 
<7c-0 The partition function characterizes the system in terms of its thermal 
properties, such as its temperature, its expectation value for energy, etc. Can 
we conclude that the physical system described by this partition function is 
a real Lorentzian black hole simply because the path integral is extremized 
by a complex (or Euclidean) black hole? The answer is yes, precisely because 
the stress-energy-momentum for the stationary complex black hole g c coin- 
cides with the stress-energy-momentum of the related stationary Lorentzian 
black hole g L . Thus, for example, when we compute the expectation value 
of energy from the partition function, it will coincide in the zero-loop ap- 
proximation with the energy surface density e of the complex black hole g c , 
which in turn characterizes the energy surface density of a real Lorentzian 
black hole g L . In this way we can conclude that the partition function indeed 
provides a description of the thermal properties of a physical black hole. 

V. Functional Integrals 

The boundary data for the partition function is chosen to be stationary on 3 B, so the 
extremum of the action can be expected to be stationary as well. 
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A path integral constructed from an action / is a functional of the quan- 
tities that are held fixed in the variational principle 51 = 0. What are held 
fixed in the variational principle are the quantities that appear varied in the 
boundary terms of 51. The fixed boundary data for the action of Eq. (4.3) 
are the metric hij on t' and t", the angle a at the corners B' and B", and 
the lapse function N, shift vector V a , and two-metric cr j, on 3 B. In the path 
integral, the gauge invariant part of the data on 3 B corresponds to the inverse 
temperature f3, chemical potential u a , and two-geometry of the boundary B. 
These are grand canonical boundary conditions. 

The inverse temperature is defined in terms of the boundary data on 3 B 
by [12] 

(3 = JdtN\ B . (5.1) 

In geometrical terms, this is the proper distance between t' and t" as mea- 
sured along the curves in 3 B that are orthogonal to the slices B. Thus, (3 
is the t-coordinate invariant part of the lapse function N. Note that (3 is a 
function on the space boundary B, as anticipated in the discussion at the 
end of Sec. 2. 

The chemical potential is defined in terms of the boundary data on 3 B 
by [12] 



Jdt V a 



Jdt N 



f dt V a \ 

B (5.2) 



JdtN\ B 



For a complex metric g c , f3u a has a geometrical interpretation as the proper 
distance along the a-coordinate line by which the spatial coordinates are 
shifted between t' and t". For a Lorentzian metric g L , u a has the following 
interpretation. Recall that the shift vector V a on 3 B is the velocity of the 
spatial coordinate system with respect to the observers who are at rest in the 
constant time slices B. It is convenient to tie the spatial coordinates to the 
local motion of the physical system [11,12], so that V a gives the local velocity 
of the system in terms of coordinate time. Thus, we see that the chemical 
potential u a is the proper velocity of the physical system as measured with 
respect to observers who are at rest at the system boundary B. In the case 
of axisymmetric boundary data, the only nonzero component of u a is the 
proper angular velocity of the system. 
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The functional integral constructed from the action / is 



p[h" , h'; a", a' 




-i[g] 



(5.3) 



where h" and h! denote the metrics on t" and t', and a" and a' denote 
the angles at the corners B" and B' . We will consider this expression as a 
functional integral over the class of complex metrics for which N, V\ and 
hij are real. 

The path integral (5.3) is the grand canonical density matrix for the 
gravitational field in a box B. The grand canonical partition function, de- 
noted Z[P,u,a], is obtained by tracing over the initial and final configura- 
tions. In the path integral language, this amounts to performing a periodic 
identification so that the manifold topology becomes M. = £ x S 1 . In ad- 
dition, a" and a' should be chosen so that the total angle a" + a! equals n. 
This insures that the boundary dM. is smooth when the initial and final hy- 
persurfaces are joined together. Thus, the grand canonical partition function 
can be written as 



The right hand side of this expression apparently depends on the angle dif- 
ference a" — a'. However, we expect that with the periodic identification, 
a" — a' is pure gauge and in a more detailed analysis would be absent from 
the path integral. 

One can consider various density matrices and partition functions corre- 
sponding to different combinations of thermodynamical variables, where one 
variable is selected from each of the conjugate pairs. For example, in ordi- 
nary statistical mechanics the thermodynamically conjugate variables might 
consist of the inverse temperature (3 and energy E, and the chemical poten- 
tial fj, and particle number N. Then the grand canonical partition function is 
Z(j3,fi), the canonical partition function is Z(/3,N), and the microcanonical 
partition function (the density of states) is v{E,N). These partition func- 
tions are related to one another by Laplace and inverse Laplace transforms, 
where each transform has the effect of switching the functional dependence 
from some thermodyanamical variable (such as (5) to its conjugate (such as 




(5.4) 



E). 
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When the gravitational field is included in the description of the system, 
all of the thermodynamical data can be expressed as boundary data [18,12]. 
In the path integral formalism, the effect of the Laplace and inverse Laplace 
transforms is simply to add or subtract certain boundary terms from the 
action. Thus, for example, the action I m appropriate for microcanonical 
boundary conditions just differs from the action / by boundary terms: 

I m = I- J 3B d 3 x^(Ne-V a 3a ) . (5.5) 

From Eq. (4.3) we see that the contribution to SI m from the boundary element 
3 B is 

SlmUs = J d 3 x (-N6(y/je) + V a 5(V^Ja) ~ {N^/2)s ab 5a ab ) . (5.6) 

The path integral constructed from I m is a functional of the quantities that 
appear varied in the boundary terms of SI m . Thus, the microcanonical den- 
sity matrix is 

p m [h",h';a",a';e,j,a] = J Dge- 1 ^ , (5.7) 

The trace of this density matrix is the density of states for the gravitational 
field in a box, 

(5.8) 

a"+a'=7r 

When combined, Eqs. (5.7) and (5.8) yield a path integral expression for 
the density of states. If the contours of integration for N and V 1 are rotated 
from the real to the imaginary axis, then iy[s,j, a] is expressed as a functional 
integral over the class of metrics with N, V 1 , and hij real. In other words, 
v[e,j, a] can be expressed as a path integral over real Lorentzian metrics [12]. 

VI. Black Hole Entropy 

In the previous section the partition function Z[f3, iv, a] and the density 
of states u[e,j,a] were constructed as functional integrals over the gravita- 
tional field on manifolds whose topologies are necessarily E x S 1 . This would 
seem to be an unavoidable consequence of deriving Z[f3,u,cr] and v[e,j,a] 



'[£,3,°] = J Dhp m [h,h;a",a';e,j,cr] 
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from traces of density matrices because the density matrices p and p m are de- 
fined in terms of functional integrals on manifolds M. with product topology 
Ex/. However, experience has shown that for a black hole, the functional 
integrals for the partition function and density of states are extremized by 
a complex (or Euclidean) metric on the manifold R 2 x S 2 [3,12]. Thus, one 
would expect the black hole contribution to the density of states to come 
from a path integral that is defined on a manifold with topology R 2 x S 2 . In 
this section we will show how the black hole density of states can be related 
to the microcanonical density matrix. 

We begin by considering the manifold M. = E x /, where E is topo- 
logically a thick spherical shell (S 2 x J). The boundary <9E = B consists of 
two disconnected surfaces, an inner sphere Bi and an outer sphere B Q . The 
boundary element 3 B consists of disconnected surfaces as well, 3 Bi = Bi x I 
and 3 B Q — B x I. The results of the previous sections can be applied in 
constructing various density matrices for the gravitational field on E. We 
wish to consider the particular density matrix that is defined through the 
path integral by the action 




i* differs from the action / of Sec. 4 by boundary terms that are not the same 
for the two disconnected parts of 3 B. The contributions to the variation 51* 
from 3 Bi and 3 B Q are 

5h\s Bi = J 3B d 3 x(( v ^e)5N-( v ^ ja )5V a + a ah 5(N v ^s ab /2)) , (6.2) 

5h\ 3Bo = [ d 3 x (-N5(yfce) + V a 5{^ 3a ) - (N^s ab /2)5a ab ) . (6.3) 

Comparison with Eq. (5.6) shows that at the outer boundary element 3 B Q 
we have chosen microcanonical boundary conditions. At the inner boundary 
element 3 Bi we have chosen completely open boundary conditions. By "com- 
pletely open" we mean that none of the traditional "conserved" quantities like 
energy, angular momentum, or area is fixed. Thus, all of the stress-energy- 
momentum quantities e, ji, and s ab are allowed to fluctuate on the inner 
boundary element while their conjugates, the inverse temperature, chemical 
potential, and spatial stress, are held fixed. The microcanonical boundary 
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conditions and the completely open boundary conditions are precisely oppo- 
site in this respect. 

Observe that the spatial stress tensor s ab can be split into its trace 
= s ab a ab = s°, which is the sum of the normal components of stress ("pres- 
sures"), and the shear stresses rf b = s ab — 9a ab /2. Then the contribution 
from 3 Bi to 51* becomes 

Sh\s Bi = J d 3 x ((V^e)5N - (V^Ja)SV a 

-a ab S(N v ^rj ab ) - ^5{N9)) . (6.4) 

Thus, the fixed data on 3 B { consist of N, V a , N^/arj ab , and N9, and the 
density matrix constructed as a path integral from is 

p,[h",h';a",a';e,j,a;N,V,N^ri,N9] = J Dge' 1 ^ . (6.5) 

Here, it is understood that the data e, j a , and are fixed on the outer 
boundary element 3 B Q and the data N, V a , y / a?7 ab , and 9 are fixed on the 
inner boundary element 3 Bi. Also, a" represents the angles at both discon- 
nected parts of B" and likewise for a'. 

The black hole density of states v*[e,j, a] is obtained from the trace of 
the density matrix along with the following special choice of data on the 
inner boundary element 3 Bf. 

N = , (6.6) 
V a = , (6.7) 
Ny/ar] ab = , (6.8) 

N9 = . y - . (6.9) 

In Eq. (6.9), (t" — t') is just the range of coordinate time t. In addition to 
the above conditions, the sum of angles a" + a! must equal ix so that the 
boundary is smooth. The resulting expression for the black hole density of 
states is 



v, 



= J dh p*[h,h; a", a'; e,j,(r;0, 0,0, -An/ n{t" -t')\ 



. (6.10) 

a"+a'=iT 
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First we will discuss the geometrical meaning of the data (6.6)-(6.9) and show 
that Eq. (6.10) indeed gives the correct expression for black hole entropy. We 
will then discuss the physical implications of this result. 

With the periodic identification, condition (6.6) fixes to zero the length 
of the circles on the inner boundary 3 Bi = B { x S* 1 that are orthogonal to the 
slices Bi. That is, N = on the inner boundary causes the hole in the planes 
defined by the unit normals and n M to become closed. This condition 
therefore seals the opening in the manifold, changing the manifold topology 
to R 2 x S 2 . The sealed inner boundary of the manifold is called the bolt [19]. 
The condition (6.7) is not mandatory, but is a convenient choice. In terms 
of a physical black hole metric, this condition fixes the spatial coordinate 
system to be co-rotating with the black hole [11,12]. The chemical potential 
is then related to the shift vector on the outer boundary element B according 
to the relationship (5.2). Without the condition (6.7), the chemical potential 
and the shift vector on B Q are not related in the simple way (5.2). Condition 
(6.8) is chosen to be consistent with condition (6.6). Finally, condition (6.9) 
ensures that the geometry is smooth at the bolt. This can be seen by writing 
Eq. (6.9) in detail. From Eq. (4.6) we have 

N 2 - 2 ^7° A-k 
- -k + -(n^N) + ^o ah -— = - , (6.11) 

which, in light of conditions (6.6) and (6.7), becomes (t" — t'){n l diN) = —2n. 
Now, (t" — t')(ri L diN) is just minus the rate of change of proper circumference 
with respect to proper radius for the circles on the inner boundary 3 Bi that 
are orthogonal to the slices B^. (The minus sign appears because n l is the 
unit normal to Bi that points outward from £.) Thus, condition (6.9) fixes 
the ratio of circumference to radius to 2ir. This insures that the spacetime 
four-geometry is smooth as the inner boundary is sealed. 

The correctness of the prescription (6.5), (6.10) can be confirmed by 
considering the evaluation of the functional integral for v*[s,j, cr], where the 
data e, j a , cr a b on the outer boundary 3 B Q correspond to a stationary black 
hole. That is, let e, j a , <J a b be the stress-energy-momentum for a topologi- 
cally spherical two-surface B a within a time slice of a stationary Lorentzian 
black hole solution g L of the Einstein equations. In the path integral for 
u^[e,j,cr], fix this data on each slice B Q of the outer boundary 3 B Q . The 
path integral can be evaluated semiclassically by searching for metrics that 
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extremize the action 7* and satisfy the conditions at both 3 B Q and 3 Bi. One 
such metric will be the complex metric g c that is obtained by substituting 
t — > —it in the Lorentzian black hole solution g L . This complex metric will 
satisfy the boundary conditions at 3 B Q because, as discussed at the end of 
Sec. 4, the substitution t — > —it does not affect the stress-energy-momentum 
quantities e, j a , and a a i } . The complex metric also will satisfy the conditions 
at 3 Bi, where B,i coincides with the intersection of the stationary time slices 
and the black hole horizon for the Lorentzian metric. This follows from the 
observation that the lapse function corresponding to the natural stationary 
time-slicing vanishes on the horizon of a stationary black hole. Likewise, 
the corresponding shift vector vanishes on the horizon for co-rotating spatial 
coordinates. Also, with an appropriate choice of period t" — t! for coordinate 
time, the complex metric will describe a smooth geometry at the bolt. (A 
single choice of period suffices for all points on B^ because n l diN is pro- 
portional to the surface gravity of the stationary black hole and is therefore 
constant across £>,.) 

In the zero-loop approximation to the path integral, the density of states 
becomes 

v*[e,j,<r] ^e- h[9c] . (6.12) 

From Eqs. (4.7) and (6.1), we see that the action 7* written in canonical form 
is 

h = J d 4 x(-iP ij h ij + NH + V'Hi) 

+ d 3 x^ (Ne - V a j a + NB/2) . (6.13) 

In evaluating this action at the solution g c , the term proportional to hij 
vanishes because of stationarity and the constraints 7i and Tii vanish because 
the equations of motion are satisfied. The remaining terms in 7* can be found 
from the boundary conditions (6.6)-(6.9), which yield 

2-7T r 

h[9c] = / d 2 x^ . (6.14) 

The integral that remains is just the area of the bolt, or the area A H of the 
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black hole event horizon. Thus, in the approximation (6.12), the entropy isQ 



S[e,j,<r] = In ^[s,j,a} « — A 



(6.15) 



With k = 8tt, this is the standard result S = A H /4 for the black hole entropy. 

It might appear that the conditions (6.6)-(6.9) have been artificially 
contrived in an effort to obtain the result (6.15). On the contrary, these 
conditions are "no boundary" conditions in the sense that they are precisely 
the conditions needed to seal the opening 3 Bi in the manifold M. = S x S , 
and convert the manifold topology to R 2 x S 2 . Thus, with these conditions 
there is no inner boundary in the spacetime four-geometry. Equations (6.6)- 

(6.9) are more properly called regularity conditions, rather than boundary 
conditions. Accordingly, the path integral defined through Eqs. (6.5) and 

(6.10) is not to be viewed as a functional of any inner boundary data. (We 
took a slightly different viewpoint in Refs. [11,12] which will be discussed 
in detail in [15].) This is why the argument of the density of states (6.10) 
includes only the conditions on the outer boundary 3 B a . 

For a general choice of data on 3 Bi, the object defined by the trace of 
the density matrix (6.5) would have to be treated as a hybrid partition func- 
tion with completely open boundary conditions on 3 B^ and microcanonical 
boundary conditions on 3 B . That is, for a general choice of data N, V a , 
^favj ah ', 9 on 3 Bi, the trace of the density matrix yields 



In order to compute the density of states in this case, we would first perform 
a series of inverse Laplace transforms to change the data on 3 Bi from com- 
pletely open conditions to microcanonical conditions. The inverse Laplace 
transforms would have the effect of canceling the boundary terms in the 
action (6.13) and changing the action to the microcanonical action I m of 
Eq. (5.5). This action is zero for any stationary solution, including station- 
ary black holes. So in the zero-loop approximation, the entropy so computed 

5 The corresponding result for an arbitrary diffeomorphism invariant theory of gravita- 
tional and matter fields has been derived in Ref. [20]. 



Z*[e,j,a;N,V a ,Vtr) ab ,8] 

= Jdhp, [h, h; a", a'; e, j, a; N, V a , ^if\ 9] 



a"+a'=Tr 



(6.16) 
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would be zero. This holds true even if the action is extremized by a black 
hole solution. 

The difference between the density of states as derived from the partition 
function Z* and the black hole density of states is precisely in the presence 
or absence of an inner boundary 3 Bi and its associated data. Thus, one is 
led to the view that the black hole entropy derived from is associated with 
a lack or absence of boundary information. The entropy derived from is 
greater than the entropy derived from because in the case of there is 
less spacetime boundary data information: less information, more entropy. 

We should remark that our present interpretation does not support the 
view [5] that the entropy derived from z/* arises from an integration over 
interior boundary data. In fact, by integrating over interior data in the den- 
sity of states derived from Z*, one merely recovers the partition function Z*. 
The black hole density of states should be viewed as the result of a path 
integral over certain gravitational fields (spacetimes) with no inner bound- 
aries whatsoever, not as the result of an integration over inner boundary 
conditions. 

The above conclusions will be developed in greater detail in forthcoming 
publications [15]. 
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